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ABSTRACT 

Future weak lensing surveys will directly probe the clustering of dark matter, 
in addition to providing a test for various cosmological models. Recent studies 
have provided us with the tools which can be used to construct the complete 
probability distribution function for convergence fields. It is also possible to 
construct the bias associated with the hot-spots in convergence maps. These 
techniques can be used in both the quasi-linear and the highly nonlinear regimes 
using various well developed numerical methods. We use these results here to 
study the weak lensing statistics of cosmological models with dark energy. We 
study how well various classes of dark energy models can be distinguished from 
models with a cosmological constant. We find that the ratio of the square root 
of the variance of convergence is complementary to the convergence skewness S3 
in probing dark energy equation of state; it can be used to predict the expected 
difference in weak lensing statistics between various dark energy models, and for 
choosing optimized smoothing angles to constrain a given class of dark energy 
models. Our results should be useful for probing dark energy using future weak 
lensing data with high statistics from galaxy weak lensing surveys and supernova 
pencil beam surveys. 

Subject headings: Cosmology: theory - weak lensing - Methods: analytical - 
Methods: statistical -Methods: numerical 



1. Introduction 

Recent cosmological observations favor an accelerating universe (Garnavich et al. 1998a; 
Riess et al. 1998; Perlmutter et al. 1999). This implies the existence of energy of unknown 
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nature (dark energy), which has negative pressure. Current data are consistent with dark en- 
ergy being a non-zero cosmological constant (see for example, Wang & Garnavich 2001; Bean 
& Melchiorri 2002). Many other alternative dark energy candidates have been considered, 
and are consistent with data as well. For example, quintessence, k-essence, spintessence, etc. 
(Peebles & Ratra 1988; Frieman et al. 1995; Caldwell, Dave, & Steinhardt 1998; Garnavich 
et al. 1998b; White 1998; Efstathiou 1999; Steinhardt, Wang, & Zlatev 1999; Podariu & 
Ratra 2000; Sahni & Wang 2000; Sahni & Starobinsky 2000; Saini et al. 2000; Waga & 
Frieman 2000; Huterer & Turner 2001; Ng & Wiltshire 2001; Podariu, Nugent, & Ratra 
2001; Weller & Albrecht 2001) 

Different dark energy models can be conveniently classified according to the equation of 
state of the dark energy component, wx- For example, for quintessence models, dwx/dz > 
0, while for k-essence models, dwx/dz < 0. There are many complimentary probes of 
dark energy. These include, the distance-redshift relations of cosmological standard candles, 
Cosmic Microwave Background Anisotropy, volume-redshift test using galaxy counts, the 
evolution of galaxy clustering, weak lensing, etc. These different methods to probe dark 
energy are complimentary to each other, and can provide important consistency checks, due 
to the different sources of systematics in each method (for example, see Kujat et al. 2002). 

Weak lensing surveys (Bacon, Refregier & Ellis 2000; Van Waerbeke et al. 2000; 
Wittman et al. 2000; Maoli et al. 2001; Van Waerbeke et al. 2001; Wilson, Kaiser, & 
Luppino 2001; Bacon et al. 2002; Hoekstra et al. 2002; Refregier, Rhodes & Groth 2002), 
currently underway and more proposed in the near future, are well suited to studying the 
dark energy equation of state. Weak lensing directly probes the gravitational clustering and 
the background cosmology. Many recent studies, both theoretical and numerical, have ana- 
lyzed these possibilities. Observational teams have already reported first detections of cosmic 
shear. On theoretical front progress has been made in modeling the statistics using both 
perturbative calculations which are valid for large smoothing angles and also using the well 
motivated hierarchical ansatz which is valid for small smoothing angles. Numerical studies 
carried out so far uses ray shooting experiments and are quite useful in testing analytical 
calculations. 

In an earlier study on probing quintessence using weak lensing, Hui (1999) concluded 
that the large scale convergence skewness can directly provide a constraint for wx, the equa- 
tion of state of dark energy. Similarly, Huterer studied the use of weak lensing convergence 
power spectrum and three-point statistics to constrain dark energy models. Both of these 
papers are useful in utilizing large weak lensing surveys of galaxies to constrain dark energy. 

In this paper we construct the complete probability distribution function of convergence 
to study the effects of dark energy. Our results apply to the weak lensing of both galaxies 
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(on large angular scales) and type la supernovae (on small angular scales). We compare two 
classes of dark energy models (one with effective constant equation of state wx, the other 
with time- varying wx) against that of a A dominated universe. 

The technique we use in this paper has been tested in detail using N-body calculations. 
Analytical results are obtained for large smoothing angles using the perturbative calculations, 
and for small smoothing angles using the hierarchical ansatz. We focus on both the one- 
point probability distribution function and the bias associated with convergence maps in 
the quasi-linear and the highly nonlinear regimes. Our studies are quite complementary 
to the studies done using a Fisher matrix analysis for the recovery of power spectra from 
observations. While such studies are well suited for recovering the nonlinear matter power 
spectra; the study of the probability distribution function and the bias will give us a direct 
probe of non-Gaussianity developed through gravitational clustering. 

The paper is organized as follows. In Section 2 we give the definition of some basic 
equations for reference, and identify the dark energy models studied in this paper. Section 
3 discusses the weak lensing statistics of dark energy models. Section 4 contains discussions 
and a summary. 



2. Notation 

The weak lensing convergence, k, maps the distribution of projected density fields, and its 
statistics is directly related to that of the underlying matter distribution. We write 

«(0o)= I*' dxw(x)8(r(x)8o,x), (1) 
Jo 

where r(x) is the angular diameter distance and iv(x) is the weight function associated with 
the source distribution. The observer is placed at x — 0, and the sources (which for simplicity 
we assume are all at the same redshift) are placed at % s . x is given by 

x(z) = cH^ 1 f dz'[n m (i + z'f + n k (i + z'f + n x /(^)]- 1/2 (2) 

Jo 

where the Q's denote the fraction of critical density in various components, fix denotes 
the dark energy component, and = 1 — Q m — Q x , The function f(z) parametrizes the 
time-dependence of the dark energy density, and f{z = 0) = 1. For dark energy with 
constant equation of state, wx = Px/px = constant, f(z) = (1 + Z )^ 1+Wx ^ . The limiting 
case with w x = —1 [he., f(z) = 1] corresponds to a cosmological constant. Note that to 
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obtain accelerated expansion of the universe, we need p + 3p < 0, which implies wx < — 1/3 
for a power law dark energy density f(z). In general, the dark energy equation of state, 
Wx(z), can be written in terms of the dimensionless dark energy density, f(z), Wx(z) = 
|(1 + z) f'(z)/f(z) — 1, for dark energy density or equation of state with arbitrary time 
dependence. (Wang & Garnavich 2001) 

We study two classes of dark energy models. The first class contains dark energy models 
with effective constant equation of state, wx = — 1/3, —2/3, —1 (ACDM), and —1.9. If dark 
energy arises from classical fields, it must satisfy the weak energy condition, which requires 
that wx > —1- The weak energy condition violating toy model, wx = —1.9, is motivated 
by quantum gravity models of inflation in which quantum effects lead to the violation of 
the weak energy condition (Onemli & Woodard 2002). The second class contains two dark 
energy toy models with linear time-varying equations of state, wx = w q (z) = — 1 + z and 
— 1 + 2z/3. These are motivated by quintessence models which have wx that effectively 
increases with z. All the models have Q m = 0.3, fix = 0.7, h = 0.7, n s = 1 (power law 
index of the primordial power spectrum), and <jg = 0.8. We normalize the nonlinear power 
spectrum to <r 8 . Table 1 lists these dark energy models. 



Table 1 

Dark Energy Models 



Model 




w x (z) 


ACDM 


-1 




constant wx 


-1/3, 


-2/3, -1.9 


w g (z) 


-l + 2z/3, 


-1 + 2 



The two classes of dark energy models are compared with the fiducial ACDM model to 
quantify the variation in various statistics of convergence maps. 

3. Statistics of Weak Lensing in Dark Energy Models 

To compute the statistics of weak lensing convergence field, we need to relate it to 
the statistics of three dimensional density field of underlying matter distribution. In recent 
studies such analysis has been done extensively. We use such a formalism to explore the 
weak lensing statistics for the dark energy cosmologies. For large smoothing angles we use 
the perturbative calculations and for small smoothing angles we use the well motivated 
hierarchical ansatz to compute the relevant quantities. 



5 



3.1. Evolution of the Matter Power Spectrum 



We compute the matter power spectrum using the scaling ansatz of Hamilton et al. 
(1991), which was later extended by various authors [see e.g. Peacock & Dodds (1996)]. 
This ansatz essentially consists of postulating that 4irk 3 P(k) = f[4irkfPi(ki)], where P(k) 
is the nonlinear power spectrum and Pi is the linear power spectrum, and the function / in 
general will depend on the initial power spectrum. The linear power spectrum is evaluated 
at a different wave number, ki — (1 + ATrk 3 P(k))~ 1 ^ 3 k, hence the mapping is non-local in 
nature. The cosmological model enters through the linear growth function g(z), so that 
Pi(k, z) = [g(z)/(l + z)] 2 Pi(k, z = 0). The linear growth function is given by: 



We compute the linear growth function by direct integration [for more on power spectrum 
evolution in quintessence models see (Benabed & Bernardeau 2001)]. 

Our method enforces stable clustering in the nonlinear regime and assumes the hier- 
archical ansatz (which is tested by numerous numerical simulations), therefore we are able 
to predict the higher order correlation functions (Davis & Peebles 1977, Groth & Peebles 
1977, Fry & Peebles 1978) and (their Fourier transforms or) the multi-spectrum correctly. 
Combining these with the powerful technique of the generating function we can construct the 
complete probability distribution function and the bias associated with convergence maps. 



Perturbative calculations depend on the expansion of the convergence field k{6q) for 
smoothing angle 9 in terms of perturbative expansion of the density field 5. Such an analysis 
can be performed in an order by order manner, 




(3) 



Where 



e{z) = y/n m (i + zf + n k (i + zf + n x /(*) 



(4) 



3.2. Convergence Probability Distribution Function 



K 




(5) 
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where 5® and «W correspond to the i-th order perturbative expansion, i — 1 being the 
linear order. In the perturbative regime at tree level (Fry 1984; Bernardeau 1992, 1994; 
Bernardeau & Schaeffer 1992), it is possible to introduce vertex generating function G{r) 
which will encode the tree level contribution from all orders. The smoothing using a top- 
hat filter function can be incorporated in the generating function formalism and then the 
generating function can be written in terms of the generating function of the unsmoothed 
case. All statistical quantities including probability distribution functions can be constructed 
once we have solved for the tree-level generating functions (Bernardeau 1992, 1994). We write 



G{t) 



Gs (r) = G 



PT 



T 



a{R,{l+G PT {r))^) 



(6) 



The parameter K a can be determined from the dynamical equation governing the evolution 
of perturbations in the quasi-linear regime, and it is given by n a = ^±§=±. The variance at a 



length scale <t(Rq) = R 



(n+3)/2 



where a local power law spectrum index n is used to evaluate 



the generating function. The generating function is now used to compute the probability 
distribution function at a particular smoothing scale. 

In the highly nonlinear regime (Balian & Schaeffer 1989, Davis & Peebles 1977, Groth 
& Peebels 1977, Szapudi & Szalay 1993, Scoccimarro & Frieman 1999, Munshi et al. 1999c), 
the perturbative series starts to diverge and the usual perturbative calculations are replaced 
by the hierarchical ansatz for higher order correlation functions, which can be built from 
two-point correlation functions. The amplitude of various contributions can be constructed 
from the knowledge of the generating function. It was found from analytical reasoning 
and numerical experimentation that the generating function in the highly nonlinear regime 
retains exactly the same form as in the quasi-linear regime; however, the value of K a is 
changed (Beranrdeau 1992) - it is now treated as a free parameter (Colombi et al. 1995, 
Munshi et al. 1999a,b). It is customary to use a different parameter uo that is easy to 
evaluate from numerical simulations, n a = ^r^y- It was found that u = .3 reproduces 
various statistics in the nonlinear regime quite well (Colombi et al. 1997, Colombi, Bouchet, 
Schaeffer 1995, Munshi et al. 1999d). To compute the probability distribution function, one 
has to compute the void probability function (f>(y), which acts as a generating function for 
normalized cumulants or Sn parameters and can be expressed in terms of the function G(t) 
as (Balian & Schaeffer 1989): 

m = yG(r) - \yr^G{r) 



(7) 
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Finally the probability distribution function can now be written as (Balian & Schaeffer 1989): 



pm r d y 



{l + 6)y-<j>(yy 



(8) 



In recent studies it was found that in both the quasi-linear and the highly nonlinear 
regimes, it is possible to introduce a reduced convergence field (Munshi & Jain 2000, Munshi 
2002), ri = K J Km '" , which to a very good approximation follows the same statistics as 1 + 8. In 
the quasi-linear regime it follows the smoothed projected density, and in the highly nonlinear 
regime it simply follows the 3D statistics of the density field. 1 



The variance of rj is given by (Valageas 2000a; Valageas 2000b) 



(9) 



with 



dxw(x,x.s), W=*J Ti~~ w2{m9o) ' (10) 

where "sinn" is defined as sinh if fi*. > 0, and sin if < 0. If = 0, the sinn and 
fife's disappear. A 2 (k,z) = 4irk 3 P(k, z), k is the wavenumber, P(k,z) is the matter power 
spectrum, The window function W(Vk6 ) = 2J 1 (Vk9 ) / (Vk9 ) for smoothing angle 9 . Here 
Ji is the Bessel function of order 1. Note that the clustering of the dark energy field would 
lead to an increase in the transfer function on very large scales. Huterer (2002) has shown 
that the clustering of the dark energy field can be neglected on the scales relevant to weak 
lensing surveys. 



Fig.l shows — K min and a/^, and Figs. 2 & 3 show the pdf for the two classes of dark 
energy models listed in Table 1. 



1 It was recently shown by Munshi (2002) that analytical results obtained by direct perturbative calcula- 
tions can also be obtained by using a functional fit obtained from assuming a log-normal evolution of local 
correlated density field. This method was also found not only to work for one point probability distribution 
function but also for bias associated with convergence maps. 
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3.3. Bias Associated with Convergence Maps 



Assuming a correlation function hierarchy guarantees that we have a two-point proba- 
bility distribution function P(ki, k 2 ) which can be factorized as follows (Munshi 2001), 



P(/Ci, K 2 )dK 1 dK 2 = P(«i)P(« 2 ) [1 + 6(ki)^ 2 6(«2)] dKidK 2 - (11) 

The function * s the two point correlation function corresponding to convergence maps. 
The function 6(/c) is the bias associated with the convergence maps, and can be shown to be 
related to the bias associated with overdense regions, b(l + 5) (Munshi 2001). The perturba- 
tive calculations also produce similar results, although the nature of the bias function 
changes from one regime to another. As was the case for one point normalized moments, 
whose generating function was related to the one point probability distribution function, the 
generating functions for two-point collapsed higher order correlation functions, also known 
as cumulant correlators, are also related to the bias function in a very similar manner. We 
write 



P°° dv 

J -too 2m 



;i + <%-</>(</) 



(12) 



where r is a generating function for the cumulant correlators. However, it turns out that 
the differential bias, as we have written down above, is difficult to estimate from numerical 
data. Therefore we work with the cumulative bias, which is the bias associated with points 
where convergence maps cross a particular threshold (Munshi 2001). Previous studies against 
numerical simulations showed that the bias function describes the numerical results quite 
accurately. It was shown in earlier studies that 6(/c) = b ( 1+s ^ . 

Figures 4 & 5 show the cumulative bias for the two classes of dark energy models 
indicated in Figure 1. Clearly, the cumulative bias of convergence is complementary to the 
convergence pdf in probing the non-Gaussianity of gravitational clustering and constraining 
dark energy models. 



3.4. A New Indicator for Deviations from the ACDM model 

We find that the deviations of dark energy models from the fiducial ACDM model can 
be quantified with a single parameter 
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l + e = 



V^(XCDM) 
vC(ACDM) 



K min (XCDM) ^(XCDM) 



« min (ACDM) ^(ACDM) ' 



(13) 



where XCDM represents an arbitrary dark energy model. Figure 6 shows the indicator 1 + e 
for the two classes of dark energy models studied in this paper, for smoothing angle 9 = V 



Comparison of Fig.6 and Figs. 2-3 shows that the more the pdf of the dark energy model 
differs from that of the fiducial ACDM model, the more the indicator 1 + e deviates from 
one. This indicates that the pdf is primarily determined by its variance, which is consistent 
with the finding of the existence of a universal probability distribution function (in terms of 
the scaled convergence if) for weak lensing amplification by Wang, Holz, & Munshi (2002). 

It is useful to compare our new indicator with the convergence skewness S3 (Hui 1999) 
for the same models. Figure 7 shows S3 for the same models as in Figure 6, with the 
same line and arrow types. We have computed the convergence skewness S3 using Hyper- 
Extended Perturbation theory in the nonlinear regime (small smoothing angular scales), 
and perturbative results are adopted for the quasi-linear regime (larger smoothing angular 
scales). While S3, as an indicator, mainly encodes the information about non-Gaussianity, 
the indicator we have proposed is directly related to the variance and is more sensitive to 
the projected density power spectrum. 

Our new weak lensing pdf shape indicator, 1 + e, is complementary to S3 in constraining 
the dark energy equation of state. The indicator, 1 + e, is sensitive to smoothing angle 9q, 
while S3 is not very sensitive to 9 at small angular scales. Feasible future supernova surveys 
can yield a large number of type la supernovae out to redshift z = 1 and beyond (Wang 
2000, SNAP 2 ). It may be possible to directly measure the weak lensing pdf with sufficiently 
high statistics (Metcalf & Silk 1999; Seljak & Holz 1999); this would allow us to utilize the 
pdf with different smoothing angles to probe different ranges of constant Wx models, and 
the variation of wx with z. 



We have analyzed weak lensing statistics for two classes of dark energy cosmological 
models. One class of dark energy models have effective constant equation of state wx, while 
the other have linear time- varying Wx(z) inspired by quintessence models. The weak lensing 



and 15'. 



4. Discussions and Summary 



2 See http://snap.lbl.gov. 
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statistics of these dark energy models are compared with that of a ACDM model. 

It has been shown that in directly using the distance-redshift relations of type la su- 
pernovae to probe dark energy, it is optimal to measure the dark energy density, px{z) — 
Px(0) f(z), instead of the dark energy equation of state, wx(z). (Wang & Garnavich 2001; 
Wang & Lovelace 2001; Tegmark 2001) In this paper, for convenience and illustration, we 
have used wx to classify various models. 

Note that we have considered a dark energy toy model which violates the weak energy 
condition, as similar models could arise from quantum effects in quantum gravity models 
of inflation (Onemli & Woodard 2002). Also, we have only considered dark energy models 
with linear time- varying wx(z), although dark energy models with much more complicated 
time dependence in wx(z) have been proposed (see for example, Bassett et al. 2002). This 
is because it is extremely difficult to extract the time dependence of wx(z), even if it were 
a simple linear function of the redshift z, from observational data (for example, see Maor et 



al. 2001; Barger & Marfatia 2001; Wang & Garnavich 2001; Kujat et al. 2002; Maor et al. 



We have studied the statistics of the cosmic convergence field via various diagnostics 
including the one point probability distribution functions and the bias associated with conver- 
gence "hot spots" . The analysis was done for both the quasi-linear scales where perturbative 
calculations are valid and also for very small angular scales where hierarchical ansatz is gen- 
erally used to quantify the statistical distributions. Following earlier studies we introduce a 
quantity K min which can help us to write the observed convergence field in terms of a reduced 
convergence field which in turn represents directly the statistics of density distribution. For 
large smoothing angles, perturbation theory predicts this quantity i] to trace the projected 
density field, and on small angular scales it traces the nonlinear density field in three di- 
mensions. The lower order moments corresponding to various cosmologies have already been 
investigated in detail and our studies complement these results. Also we have used top-hat 
window functions, but Bernardeau & Valageas (2000) have shown how to generalize similar 
calculations in terms of aperture mass using compensated filters. 



which can be used to predict the expected difference in weak lensing statistics between various 
dark energy models, and for choosing optimized smoothing angles to constrain a given class 
of dark energy models. For example, small smoothing angles are favored for constraining 
dark energy models with wx < — 1. 

Our proposed 1 + e statistics is related to the volume average of two-point statistics 
of £ K . Note that while t; v only depends on the smoothing angle and the underlying mass 



2002). 



We have identified a new weak lensing pdf shape indicator, 1 + e oc 
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distribution, n min encodes the dependence on cosmological parameters. Therefore, our new 
statistical indicator, 1 + e oc = K min y/Cqj is of interest as far as we are interested in 
differentiating various cosmological models. We found that both P(k) and b(K) depends 
on k m i n and £ T; , however, it is difficult to infer the difference in these statistics in various 
cosmologies. The 1 + e statistics we have devised, however, is quite interesting in the sense 
that it can very easily be used to check how much various dark energy models differ in weak 
lensing. It can be used to complement and supplement various other statistics such as S3. 

We have computed the convergence skewness S3 for the dark energy models considered 
in this paper. We find that the new weak lensing pdf shape indicator, 1 + e, is indeed 
complementary to S3 in probing dark energy equation of state. 

We note that getting maps of convergence is difficult compared to direct evaluation of 
non-Gaussian statistics from shear maps (see e.g. Schneider & Lombardi 2002; Zaldariaga 
& Scoccimarro 2002). On the other hand, the construction of convergence statistics can be 
directly modeled at arbitrary level, whereas for shear field the computation of statistics is 
done in a order by order manner so far. So an independent analysis of convergence maps 
constructed from shear maps should be useful in constraining various errors which might get 
introduced during various stages of data reduction. The question of error bars in weak lensing 
measurements has been dealt with in great detail in Munshi & Coles (astro-ph/0003481, to 
appear in MNRAS) for various window functions and are independent of the cosmological 
model assumed. Our convergence statistics can be a powerful diagnostic and complementary 
tool to the shear map statistic. 

With high statistics data from future weak lensing surveys of galaxies and supernova 
pencil beam surveys (Wang 2000; SNAP), weak lensing can be a useful tool in differentiating 
different dark energy models. 

DM was supported by PPARC grant RG28936, and YW was supported in part by NSF 
CAREER grant AST-0094335. DM would like to thank Alexandre Refregerier for many 
useful discussion, and Francis Bernardeau for making a copy of his code to compute the pdf 
and bias available to us. 
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Fig. 1. — The minimum value of the convergence field, K m i n , and variance of the scaled 
convergence r\ = 1 + K,/\K min \, (, v , as functions of constant dark energy equation of state, 
wx, for source redshift z s = 0.5 and 1.0. The arrows indicate the —n m in and ^J~£^ values for 
dark energy models with time-varying equation of state, w q (z) = — 1 + z (long arrows) and 
Wg(z) = —l + 2z/3 (short arrows). Note that K min does not depend on the smoothing angle 
6q but it depends on the background dynamics of the universe. 

Fig. 2. — PDF associated with constant wx models. Two different regimes are considered 
for computing the PDF. Upper panels 9 S = l',4' correspond to the nonlinear calculations 
where we have assumed a hierarchical ansatz for the correlation hierarchy. For the lower 
panels perturbative results are used to construct the PDF at smoothing angles 9 S = 15' 
and 60'. Various curves correspond to various quintessence models. In each panel solid 
line represents the ACDM model, short dashed line represent wx = —2/3, dot-dashed line 
represents wx = —1/3, and long dashed line represents wx = —1.9. 

Fig. 3. — PDF associated with time- varying wx models compared to ACDM. As in the 
previous figure two different regimes are considered for computing the PDF. Upper panels 
9 S = 1', 4' correspond to the nonlinear calculations where we have assumed a hierarchical 
ansatz for the correlation hierarchy. For the lower panels perturbative results are used to 
construct the PDF at smoothing angles 9 S = 15' and 60'. Various curves correspond to 
various quintessence models. In each panel solid line represents the ACDM model, short 
dashed line represent wx = —I + z model, and long dashed line represents wx = — 1 + 2z/3 
model. 

Fig. 4. — Bias associated with constant w x models. Two different regimes are considered 
for computing the PDF. Upper panels 9 S = l',4' correspond to the nonlinear calculations 
where we have assumed a hierarchical ansatz for the correlation hierarchy. For the lower 
panels perturbative results are used to construct the PDF at smoothing angles 9 S = 15' 
and 60'. Various curves correspond to various quintessence models. In each panel solid 
line represents the ACDM model, short dashed line represent wx = —2/3, dot-dashed line 
represents wx = —1/3, and long dashed line represents wx = —1.9. 

Fig. 5. — Bias associated with time-varying wx models compared to ACDM. As in previous 
figure two different regimes are considered for computing the PDF. Upper panels 9 S = l',4' 
correspond to the nonlinear calculations where we have assumed a hierarchical ansatz for 
the correlation hierarchy. For the lower panels perturbative results are used to construct the 
PDF at smoothing angles 9 S = 15' and 60'. Various curves correspond to various quintessence 
models. In each panel solid line represents the ACDM model, short dashed line represent 
w x = —1 + z model, and long dashed line represents w x = —1 + 2z/3 model. 
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Fig. 6. — Indicator of the deviations of dark energy models from the fiducial ACDM model, 
1 + e, as function of constant dark energy equation of state wx, for smoothing angle 6q = 1' 
and 15'. The arrows indicate the 1 + e values for dark energy models with time-varying 
equation of state, w q (z) = — 1 + z (long arrows), and w q (z) = — 1 + 2z/3 (short arrows). 

Fig. 7. — The skewness S3, for the same models as in Fig.6. The line and arrow types are 
the same. In the left panel we use Hyper-Extended Perturbation theory to compute the 
convergence skewness, whereas in the right panel (where larger smoothing angular scales are 
considered) perturbative results are adopted. 
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Fig. 1. — The minimum value of the convergence field, n m i n , and variance of the scaled 
convergence i] = 1 + K>/\K m in\, Cr?> as functions of constant dark energy equation of state, 
Wx, for source redshift z s = 0.5 and 1.0. The arrows indicate the —n m in and values for 
dark energy models with time-varying equation of state, w q (z) = — 1 + z (long arrows) and 
Wg(z) = —l + 2z/3 (short arrows). Note that K min does not depend on the smoothing angle 
6*0 but it depends on the background dynamics of the universe. 
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Fig. 2. — PDF associated with constant wx models. Two different regimes are considered 
for computing the PDF. Upper panels 9 S = l',4' correspond to the nonlinear calculations 
where we have assumed a hierarchical ansatz for the correlation hierarchy. For the lower 
panels perturbative results are used to construct the PDF at smoothing angles 6 S = 15' 
and 60'. Various curves correspond to various quintessence models. In each panel solid 
line represents the ACDM model, short dashed line represent w x = —2/3, dot-dashed line 
represents wx = —1/3, and long dashed line represents wx = —1.9. 
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Fig. 3. — PDF associated with time- varying wx models compared to ACDM. As in the 
previous figure two different regimes are considered for computing the PDF. Upper panels 
9 S = 1', 4' correspond to the nonlinear calculations where we have assumed a hierarchical 
ansatz for the correlation hierarchy. For the lower panels perturbative results are used to 
construct the PDF at smoothing angles 9 S = 15' and 60'. Various curves correspond to 
various quintessence models. In each panel solid line represents the ACDM model, short 
dashed line represent wx = —1 + z model, and long dashed line represents wx = —1 + 2z/3 



model. 




Fig. 4. — Bias associated with constant wx models. Two different regimes are considered 
for computing the PDF. Upper panels 9 S = 1', 4' correspond to the nonlinear calculations 
where we have assumed a hierarchical ansatz for the correlation hierarchy. For the lower 
panels perturbative results are used to construct the PDF at smoothing angles 6 S = 15' 
and 60'. Various curves correspond to various quintessence models. In each panel solid 
line represents the ACDM model, short dashed line represent wx = —2/3, dot-dashed line 
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represents wx = —1/3, and long dashed line represents wx = —1.9. 
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Fig. 5. — Bias associated with time-varying w x models compared to ACDM. As in previous 
figure two different regimes are considered for computing the PDF. Upper panels 9 S = 1', 4' 
correspond to the nonlinear calculations where we have assumed a hierarchical ansatz for 
the correlation hierarchy. For the lower panels perturbative results are used to construct the 
PDF at smoothing angles 8 S = 15' and 60'. Various curves correspond to various quintessence 
models. In each panel solid line represents the ACDM model, short dashed line represent 
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wx = —1 + z model, and long dashed line represents wx = — 1 + 2z/3 model. 




Fig. 6. — Indicator of the deviations of dark energy models from the fiducial ACDM model, 
1 + e, as function of constant dark energy equation of state wx, for smoothing angle #o = 1' 
and 15'. The arrows indicate the 1 + e values for dark energy models with time-varying 
equation of state, w q (z) = — 1 + z (long arrows), and w q (z) = — 1 + 2z/3 (short arrows). 
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Fig. 7. — The skewness S3, for the same models as in Fig.6. The line and arrow types are 
the same. In the left panel we use Hyper-Extended Perturbation theory to compute the 
convergence skewness, whereas in the right panel (where larger smoothing angular scales are 
considered) perturbative results are adopted. 



